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Abstract 

The notion of fluctuation indices, characterizing thermodynamic stability of sta- 
tistical systems, is advanced. These indices are especially useful for investigating the 
stability of nonuniform and trapped atomic assemblies. The fluctuation indices are cal- 
culated for several systems with Bose-Einstein condensate. It is shown that: the ideal 
uniform Bose-condensed gas is thermodynamically unstable; trapped ideal gases are 
stable for the confining dimension larger than two; trapped gases, under the confining 
dimension two, are weakly unstable; harmonically trapped gas is stable only for the 
spatial dimension three; one-dimensional harmonically trapped gas is unstable; two- 
dimensional gas in a harmonic trap represents a marginal case, being weakly unstable; 
interacting nonuniform three-dimensional Bose-condensed gas is stable. There are no 
thermodynamically anomalous particle fluctuations in stable Bose-condensed systems. 
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1 Introduction 



Thermodynamic stability of statistical systems is an important notion characterizing the 
possibility of existence of equilibrium systems as such. There are several stability conditions 
that are required to be fulfilled in order that the system be thermodynamically stable [1]. 
In the present paper, we concentrate on the stability related to fluctuations of observable 
quantities. These fluctuations are to be thermodynamically normal as soon as one assumes 
that the considered system is in thermal equilibrium. In the other case, if the fluctuations 
of at least one of the observables are thermodynamically anomalous, this implies that this 
observable cannot be measured and the system equilibrium is actually destroyed by such 
fluctuations, which explains their naming as anomalous. 

The problem of instability, caused by fluctuations, has recently attracted great attention 
with respect to nonuniform confined systems, such as trapped atomic gases. Especially 
intensive discussions on the type of particle fluctuations have accompanied the study of 
systems with Bose-Einstein condensate. Description of the main properties of the latter 
systems can be found in the book [2] and review articles [3-10]. 

In discussions of particle fluctuations in Bose-condensed systems, there have been the 
widely spread heresy that these fluctuations could be anomalous, been drastically different 
from such fluctuations in other systems. This controversy has been investigated in detail 
in Refs. [4,10-13], where it has been explained that the appearance of such anomalous 
fluctuations is merely due to calculational errors. 

When analyzing whether fluctuations of observable quantities are thermodynamically 
normal or anomalous, it is necessary to resort to the notion of thermodynamic limit. In 
the present paper, a novel characteristic is introduced, allowing for a convenient quantitative 
description of the type of the limiting behaviour of fluctuations under the increasing number 
of particles. This is the fluctuation index associated with the considered observable. The 
value of this index defines when the fluctuations are thermodynamically normal and when are 
not, thus, showing whether the system is thermodynamically stable or unstable. The notion 
is illustrated by calculating the fluctuation indices for different Bose-condensed systems: 
confined in a box, trapped by power-law potentials, and generally nonuniform. 



2 Fluctuation indices 

2.1 Definition and stability conditions 

The standard definition of thermodynamic limit implies that the number N of particles in 
the system and its volume V tend to infinity so that 

N 

N — > oo , V — > oo , — — > const . (1) 

This definition assumes that there is a well defined volume of the system. 

However, for systems confined in trapping potentials, the volume may be not fixed. In 
such a case, it is necessary to resort to a more general definition of thermodynamic limit. 
The latter can be done as follows [10,14]. Let An be an extensive observable quantity for 
a system with the number of particles N. Then the thermodynamic limit, related to this 
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observable, is defined as 



An 

N — > oo , A N — > oo , — — > const . (2) 
From here, it is seen that 

In I Am I 

lim 1 Nl < 1 . 3 

In the particular case, when the system volume V is well defined, hence the extensive ob- 
servable A^ oc V, definition (2) reduces to the standard form (1). 

Fluctuations, associated with an observable, represented by a self-adjoint operator A, are 
described by the dispersion, or variance, of this operator, given by the statistical average 

var(i) = ((i-(i)) 2 ) = (i 2 )-(i) 2 . 

The fluctuation index, related to the observable quantity (A), represented by a self-adjoint 
operator A, is defined as the limit 

i 2\ i- lnvar(A) 
(pi A) = hm — . (4) 

Since the variance itself is an extensive observable, it has to satisfy the stability conditions 

, In , var(A) . . 

<p(A)<l, hm — ^<oo, (5) 

in agreement with Eq. (3). The first of these conditions is necessary, characterizing the 
behavior of fluctuations in thermodynamic limit, while the second condition is necessary 
and sufficient. When both these conditions are valid, the system is stable. It may happen 
that the first condition is valid, but the second is not. In the latter case, we shall say that 
the system is weakly unstable. 

Fluctuations, satisfying Eq. (5), are called thermodynamically normal, while fluctuations, 
for which conditions (5) are not valid, are termed thermodynamically anomalous. This is 
because an observable with normal fluctuations can be measured, but that with anomalous 
ones, cannot, as soon as the magnitude of fluctuations, hence the measurement uncertainty, 
is larger than the observable itself. A statistical system, where at least one of the observ- 
ables exhibits anomalous fluctuations, cannot be in thermal equilibrium. Such a system is 
thermodynamically unstable. 

2.2 Uniform ideal gas 

Let us consider particle fluctuations in a system with Bose-Einstein condensate. The particle- 
number operator N is the sum 

N = N Q + jVi 

of the terms corresponding to condensed particles (No) and uncondensed particles (Ni). In 
thermodynamic limit, by the Bogolubov theorem [15], the operator Nq becomes a nonoper- 
ator number, so that 

var(iVo) ^0 (N ->■ oo) . 
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This fact is discussed in detail in review [8]. Therefore, particle fluctuations are completely 
due to uncondensed particles, 

var(iV) = var(A\). (6) 

For the ideal uniform Bose-condensed gas, confined in a volume V, the condensation 
temperature is 

2/d 

(7) 



nn 27T 

T — 

J- c — 

m 



P 



C(d/2) 



where d is space dimensionality, p is average density, and £(•) is the Riemann zeta function. 
Below this temperature, the condensate fraction behaves as 

, rp . d/2 

no = 1 - [ ¥ J (T < T c ) . (8) 

However, a finite critical temperature (7) exists only for d > 3. This is a general feature 
of the absence of spontaneous breaking of continuous symmetry in low-dimensional systems 
[16]. 

In the three-dimensional space, one has 

var (iV 1 )= (^T)V 4 / 3 . (9) 

Taking into account Eq. (6) gives the fluctuation index 

<p(N) = i (10) 

that contradicts the stability condition (5). This tells us that the ideal uniform Bose- 
condensed gas is unstable, being a pathological object with thermodynamically anomalous 
particle fluctuations. 



2.3 Trapped Bose gas 

An interesting question is whether the ideal Bose-condensed gas could be stabilized being 
trapped in an external potential. The most often used shape of the trapping potential is of 
the power-law form 

rl 

L 



U(r) = Ey 

Q=l Z 



(11) 



in which n a > and the trap frequencies and characteristic lengths are connected by the 
relations 

1 1 

= —pr , la = • (12) 

It is convenient to introduce the effective trap frequency and length by the geometric averages 

/ d x 1 /^ 1 , d x 1 /^ 1 

"°=[MH = ^i' ''"15/") (13) 
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Another important quantity, defining the confining power of potential (11), is the con- 
fining dimension 

d d 1 

•-s + £=- (14) 

where d is the real-space dimension. 

The properties of the ideal trapped Bose gas can be accurately described by the gener- 
alized quasiclassical approximation [14]. Bose-Einstein condensation occurs at the critical 
temperature 

1/8 

T c = ^jtt , (15) 



Bg s (l) 
in which 

2* JL r(l + l/n a ) 

B = 11 l/2+l/n Q ' ^ 16 ) 

and the generalized Bose function 

, s 1 f°° ZU 3 " 1 , 

9s{z) = W)L — du (17) 

is introduced [14], with the lower limit in the integral being 

f^- (18) 

Note that the standard Bose function corresponds to the limiting case of Uq — 0. Below the 
critical temperature (15), the condensate fraction is 

n = l-(£f (T<T C ). (19) 

But the formal occurrence of a critical temperature does not necessarily mean the real 
existence of a stable Bose-condensed system. To check the stability with respect to particle 
fluctuations, we have to calculate the related fluctuation index. The trapping potential 
(11) extends to infinity, so that the system volume is not fixed. Hence, the general form 
of thermodynamic limit (2) is to be employed. As an extensive quantity, we can take the 
internal energy En, considering the thermodynamic limit in the form 

N ->■ oo , E N ->■ oo , ^ ->■ const . (20) 
For the internal energy, we find 

E N = Bsg 1+S (1)T 1+S , (21) 

which transforms limit (20) into 

B , x 

N ->■ oo , B ->■ oo , — -> const . (22) 
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For the usual case of unipower trapping potentials, when n a = n, the confining power 
(14) is 

> = {t + k) d - < 23 > 

And quantity (16) becomes 

B = ^/ d (^ + -)- (24) 



7i d / 2 Un V n 



Then the thermodynamic limit (22) reduces to 

N ->■ oo , w -> , Nuq ->■ const . (25) 

Looking at limit (22), there arises a temptation to treat the quantity B as an effective 
volume. The latter, however, is not uniquely defined. And the most important is that such 
a quantity B cannot be used as a thermodynamic variable. Attempting to use it as such 
would lead to inconsistent thermodynamic relations. That is, though B reminds something 
like an effective volume, there is no any sense of identifying it with the latter. 

Under the thermodynamic limit (22), or (25), the critical temperature (15) behaves as 

T c oc A^ 1/s ^ (s < 1) , 

T c oc (lnA^)" 1 ->■ (s=l), 

T c ->■ const (s > 1) . (26) 

For the confining dimension s < 1 the critical temperature tends to zero. Hence, only s > 1 
provides a finite critical temperature. 

Again, the occurrence of a condensation temperature does not guarantee the existence of 
a stable Bose-condensed system. We need to find the fluctuation indices. For the variance 
of the particle number, we get 

var(A0 = ^#(;DV (27) 



This becomes negative for s < 1, which contradicts the definition of the variance as a non- 
negative quantity. Thus, only s > 1 can be considered. With the generalized Bose function 
(17), we find 
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var(A0 = 2 — (s = 1) , 

N /2T \ 2-5 (T \ 2 

var(A0 = - ttttw, r — \tA (1<s<2), 

1 ; (2 - s)C(s)r(s - 1) \cu J \TJ 1 ; ' 

N /T\ 2 , /2T f 



var(iV) = ^ S ~^ ( N (s > 2) . (28) 



C(«) VT C 

This yields the fluctuation indices 



= 2 (a = 1) 
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<p(N) = - (Ks<2), 
s 

<p(N) = 1 + ( S = 2), 
<p(N) = 1 (s > 2) , (29) 

where the notation 

lnlniV 
lim - Ar = +0 

n^oo In N 

is used. 

Consequently, ideal trapped gas can form an absolutely stable Bose-condensed system 
only for s > 2. The case s = 2 is on the boundary of stability. Strictly speaking, var(TV) 
diverges as N tends to infinity, but this divergence is weak, being of logarithmic type. This 
means that the Bose-condensed system, with the confining dimension s = 2, is weakly 
unstable. 

Remembering definition (14) gives the stability condition 

d d 1 

i + S^ 2 - (30) 

2.4 Harmonic trapping potential 

The most commonly considered shape of trapping potentials is that of harmonic potential, 
when n a = 2. Then s = d and B = l/u d . The condensation temperature (15) gives 

= \n{2N) {d=lh 



T c = cu 



N 



i/d 

(d > 2) . (31) 



One often states that the formal existence of the critical temperature T c implies the 
possibility of getting Bose-Einstein condensate in one- and two-dimensional harmonic traps. 
But, as has been stressed above, the mere occurrence of T c does not guarantee that such 
a Bose-condensed system would be stable, hence, could really exist. We have to check the 
system stability. 

The thermodynamic limit (25), for harmonic traps, takes the form 

N ->• oo , uj ->■ , Nuq ->■ const . (32) 

Variance (27) yields 

var(iV) = 2 (—) (d = 1) , 

wo/ 



,U J \UJ J 
7T 2 N (T 

603) \T C 



7T 2 N /T\ 3 

var(iV) = -— - (d = 3). (33) 
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As a result, the fluctuation indices are 

<P(N) = 2 (d=l), 

(p(N) = 1 + (d = 2) , 

<P(N) = 1 (d = 3) , (34) 

with the same notation for +0 as above. 

This tells us that, according to the stability condition (5), only a three-dimensional 
harmonic trap can house stable ideal Bose-condensed gas. Bose-Einstein condensation cannot 
occur in one-dimensional harmonic traps. And two-dimensional gas in a harmonic trap is 
weakly unstable. 



2.5 Interacting nonuniform gas 

Let us consider an arbitrary nonuniform Bose system of atoms interacting through repulsive 
forces. The general expression, characterizing particle fluctuations, is given [4,17] by the 
variance 

var(iV) = N+ I p(r)p(r')[g(r, r') - 1] drdr' , (35) 
where the total particle density 



p(r) = p (r) + pi(r) 



(36) 



is the sum of the condensate density po( r ) and the density of uncondensed particles pi(r), 
and g(r, r') is the pair correlation function. This expression is valid for any system whether 
equilibrium or not. 

For an equilibrium system, we shall use the local-density approximation [17] in the frame 
of the self-consistent mean-field approach [18-22]. Then the density of uncondensed particles 
can be written [17] in the form 



Pi 



(r) = fn(k,. 



in which 



n(k, r) 



u(k, r) 



2e(k,r) 

is the local momentum distribution, the notation 



dk 

(27p : 

g(k,r) 
2T 

k 2 

u(k,r) = mc 2 (r) + — 
2m 



(37) 



coth 



1 

2 



is used, and 



e(k,r) 



\ 



' k 2 V 
2m) 



c 2 (r)k 2 + 

is the local Bogolubov spectrum. The local sound velocity is defined by the equation 

mc 2 (r) = [po(r)+<7i(r)]$ , (38) 
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where the interaction strength is given by 



$ = / $(r) dr = 4tt — > , 
7 m 



with positive scattering length a s . And the anomalous average 

dk 



(r) = /^k,: 



is expressed through 



a(k, r) = 



mc 2 (r) 



coth 



(2tt) 3 
efk, r) 



(39) 



2T 



2e(k,r) 

Substituting in the right-hand side of variance (35) the Bogolubov shift [15] and keeping the 
terms up to second order with respect to the operators of uncondensed particles [4,17], we 
have 

var(iV) = N + 2 lim / p(r) [n(k, r) + <r(k, r)] dr . 
Accomplishing here the limit k — > yields 



var(iV) = T - [ ^ dr 



m J cr r 



(40) 



Bose-Einstein condensation is accompanied by the global gauge symmetry breaking, the 
latter being the necessary and sufficient condition for the former [8,23]. When Bose-Einstein 
condensate is present, the sound velocity c(r) is nonzero. Then the integral in Eq. (40) 
is proportional to N. Therefore the fluctuation index, describing particle fluctuations, is 
4>(N) = 1. The latter means that for an arbitrary nonuniform three-dimensional system of 
repulsive atoms, with Bose-Einstein condensate, particle fluctuations are thermodynamically 
normal. 



3 Conclusion 

The notion of fluctuation indices for operators, representing observable quantities, is in- 
troduced, characterizing the stability properties of statistical systems. Thermodynamically 
anomalous fluctuations imply the system instability, while stable systems exhibit thermo- 
dynamically normal fluctuations. The notion is illustrated by calculating the fluctuation 
indices for the number-of-particle operator for systems with Bose-Einstein condensate. It is 
shown that the ideal uniform Bose-condensed gas is unstable. The ideal gas can be stabi- 
lized in trapping potentials, provided that the confining dimension is larger than two. The 
case of the confining dimension two is marginal, corresponding to a weakly unstable system. 
Atomic gases in harmonic traps are stable only in the three-dimensional case. Trapped one- 
dimensional gases are unstable, and trapped two-dimensional gases are weakly unstable. A 
three-dimensional Bose-condensed system of atoms, interacting through repulsive forces, is 
stable for arbitrary external potentials. 

The particular nature of a condensate can be different. For short, we have been talking 
of atomic condensates. But there exist now several types of molecular condensates (see 
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Refs. [24-28] and review [10]). The above consideration can be applied for any type of Bose 
condensates, whether atomic or molecular. 

Only equilibrium systems have been considered in the paper. Fluctuations in nonequi- 
librium systems is a different topic (see, e.g., [29-35]). But the notion of fluctuation indices 
can be applied to nonequilibrium systems as well, since the definition of fluctuation indices 
(4) is equally valid for any type of averages, whether over equilibrium or nonequilibrium 
ensembles. The principal difference between equilibrium and nonequilibrium systems is that 
the stability conditions (5), generally, are not required for the latter. The strength of fluctu- 
ations in nonequilibrium systems can be of arbitrary magnitude. Though some restrictions 
on the fluctuation strength could be connected with steady and quasiequilibrium states. 

The fluctuation indices, introduced in the present paper, provide a quantitative charac- 
teristic of fluctuation strength, associated with the operators of observables. The knowledge 
of these indices can help for deciding under what conditions and for what kind of traps one 
could realize stable atomic systems in experiment. 
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